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Abstract. A superintegrable, discrete model of the quantum isotropic oscillator in two-dimensions 
is introduced. The system is defined on the regular, infinite-dimensional N x N lattice. It is governed 
by a Hamiltonian expressed as a seven-point difference operator involving three parameters. The 
exact solutions of the model are given in terms of the two-variable Meixner polynomials orthogonal 
with respect to the negative trinomial distribution. The constants of motion of the system are 
constructed using the raising and lowering operators for these polynomials. They are shown to 
generate an su(2) invariance algebra. The two-variable Meixner polynomials are seen to support 
irreducible representations of this algebra. In the continuum limit, where the lattice constant tends 
to zero, the standard isotropic quantum oscillator in two dimensions is recovered. The limit process 
from the two-variable Meixner polynomials to a product of two Hermite polynomials is carried out 
by involving the bivariate Charlier polynomials. 
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1. Introduction 

The purpose of this paper is to present a discrete model of the two-dimensional quantum 
oscillator that is both superintegrable and exactly solvable. The wavefunctions of 
this system will be given in terms of the two-variable Meixner polynomials and the 
constants of motion will be seen to satisfy the 5u(2) algebra. 

A considerable amount of literature can be found on superintegrable systems and 
there is a sustained interest in enlarging the documented set of models with that 
property. Recall that a quantum system with d degrees of freedom governed by a 
Hamiltonian H is said to be superintegrable if it possesses, including H itself, 2d-l 
algebraically independent constants of motion, that is operators that commute with 
the Hamiltonian. The quintessential example of a quantum superintegrable system 
is the two-dimensional harmonic oscillator, whose constants of motion generate the 
5u(2) algebra. One of the motivating observations behind the study of superintegrable 
systems is that they are exactly solvable, which makes them prime candidates for 
modeling purposes. Also of importance is the fact that these systems form a bedrock 
for the analysis of symmetries, of the associated algebraic structures and their 
representations, and of special functions. The majority of quantum superintegrable 
models cataloged so far comprises continuous systems, but there has also been some 
progress in the study of discrete systems [19, 20]; for a review on superintegrable 
systems (mostly continuous ones) and their applications, see [21]. 

In the past years, several discrete models of the one-dimensional quantum 
oscillator, either finite or infinite, were introduced [2, 5, 15, 16]. The most studied 
system, originally proposed in [6] as a model of multimodal waveguides with a finite 
number of sensor points, has 5u(2) as its dynamical algebra. In this model, the 
Hamiltonian, the position and momentum operators are expressed in terms of 5u(2) 
generators, the eigenstates of the system are the basis vectors of unitary irreducible 
representations of 5u(2) and the wavefunctions are expressed in terms of the one- 
variable Krawtchouk polynomials. As a result, the Hamiltonian has a finite number of 
eigenvalues and the spectra of the position and momentum operators are both discrete 
and finite. Germane to the present paper is also the discrete oscillator model based on 
the univariate Meixner polynomials and related to the 5u(l, 1) algebra considered in [2]. 
See also [17], where instead the Meixner-Pollaczek polynomials are involved. 

The Krawtchouk one-dimensional finite/discrete oscillator has been exploited to 
construct finite/discrete systems in two dimensions. Two approaches have been used. 
The first approach consist in taking the direct product of two one-dimensional 5u(2) 
systems to obtain a system defined on a square grid with 5u(2) ® 5u(2) as its dynamical 
algebra [3]. In the second approach [4], the isomorphism 5o(4) = 5u(2)®5u(2) is exploited 
to obtain a description of the finite oscillator on the square grid in terms of discrete 
radial and angular coordinates. In the continuum limit, both of these models tend to the 
standard two-dimensional oscillator. However, they do not exhibit the 5u(2) invariance, 
or S5rmmetry algebra, of the standard two-dimensional oscillator. 
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Recently another discrete and finite model of the two-dimensional oscillator was 
proposed [20]. This model is defined on a triangular lattice of a given size and, like the 
standard oscillator in two dimensions, it is superintegrable and has 5u(2) for symmetry 
algebra. The wavefunctions of the model, which support irreducible representations of 
5u(2) at fixed energy, are given in terms of the two-variable Krawtchouk polynomials 
introduced by Griffiths [11]. These polynomials of two discrete variables are orthogonal 
with respect to the trinomial distribution. As required, this model tends to the standard 
two-dimensional oscillator in the continuum limit. 

Shortly after [20] appeared, it was recognized that the d-variable Krawtchouk 
polynomials of Griffiths arise as matrix elements of the unitary representations of the 
rotation group SO{d + 1) on oscillator states [10]. This interpretation has provided 
a cogent framework for the characterization of these orthogonal functions and has 
led to a number of new identities. The group theoretical interpretation was also 
extended to two other families of discrete multivariate polynomials: the multivariate 
Meixner and Charlier polynomials, orthogonal with respect to the negative multinomial 
and multivariate Poisson distributions, respectively. The multi-variable Meixner 
polynomials, also introduced by Griffiths [12], were shown to arise as matrix elements 
of unitary representations of the pseudo-rotation group SO{d,l) on oscillator states 
[9]. As for the multivariate Charlier polynomials, they were first introduced as matrix 
elements of unitary representations of the Euclidean group on oscillator states [8]. Let 
us note that these family of multivariate polynomials also arise in probability theory in 
connection with the so-called Lancaster distributions [13]. 

In this paper, we present a new discrete oscillator model in two-dimensions based 
on the two-variable Meixner polynomials. The model is defined on the regular infinite¬ 
dimensional N X N lattice. It is governed by a Hamiltonian involving three independent 
parameters expressed as a 7-point difference operator. This operator is obtained by 
combining the two independent difference equations satisfied by the bivariate Meixner 
polynomials. By construction, the wavefunctions of the model are given in terms of 
these two-variable polynomials. The energies of the system are given by the non¬ 
negative integers N = 0,1,2,... and exhibit a {N + l)-fold degeneracy. Using the raising 
and lowering relations for the two-variable Meixner polynomials, the constants of 
motion of the system are constructed and are shown to close onto the 5u(2) commutation 
relations. In the continuum limit, in which the lattice parameter tends to zero, the 
model contracts to the standard quantum harmonic oscillator, as required for a discrete 
oscillator model. The contraction process is illustrated at the wavefunction level using 
the two-variable Charlier polynomials in an intermediary step. The continuum limit is 
also displayed at the level of operators. 

Here is the outline of the paper. In Section two, the essential properties of the 
two-variable Meixner polynomials are reviewed. In Section three, the Hamiltonian of 
the model is defined, the constants of motion are constructed, and the wavefunctions 
are illustrated. In Section four, the continuum limit of the model and wavefunctions is 
examined. We conclude with an outlook. 
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2. The two-variable Meixner polynomials 

We now review the properties of the two-variable Meixner polynomials using the 
formalism and notation developed in [9]. Let /3 ^ 0 be a positive real number and let 
A G 0(2, 1) be a 3 X 3 pseudo-rotation matrix. This implies that A satisfies 

T] A - jj, 

where rj - diag(l, 1,-1) and where A^ denotes the transposed matrix. In general, A 
can be parametrized by three real numbers akin to the Euler angles. The two-variable 
Meixner polynomials, denoted by are defined by the generating function 

' All Ai 2 A 21 A 22 Y^(-l '^31 A 32 

= z E (2.1) 

ni=0 112=0 V 

where the Aij are the entries of the parameter matrix A and where (j 6 )„ stands for the 
Pochhammer symbol defined as 


()6)n = 


1 n-Q 

YYllliP + k) n = 1,2,3,... 


It can be seen from (2.1) that MY^ j^ 2 ^x\,X 2 ) are polynomials of total degree n\ + n 2 in 
the variables xi and X 2 . The functions MY‘l^J^J,x\,X 2 ) satisfy the orthogonality relation 

00 00 

1 ^ 2 ) — ’ ( 2 . 2 ) 

X1 = 0X2=0 1’ 2 1 2 

where a»(xi,X 2 ) is the negative trinomial distribution 

e;(xi,X 2 ) = (1 - Cl - C2f cY4'’ (2.3) 

xi!x 2 ! 

and where the parameters ci, C 2 are given by 


The polynomials MII^^^^(xi,X 2 ) have an explicit expression in terms of Aomoto-Gelfand 
hypergeometric series [14]. This expression reads 


M|f/,„,(X1,X2) = (-I)"l^'^2^ 


(P)ni+n 2 ( Asi)"^^ ( A 32 


= E 

/J,V,p,( 7>0 


y ni!n 2 ! vA 33 j (A 33 
(—ni)^+v(—n2)p+(7(—xi)^+p(—X2 )v+(7 

plvlplal il3)^+y+p+a 


(1-Uii)^(1-U2i)^(1-U12)^(1“W22)° 


(2.4) 
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where the parameters Uij are given by 

A 11 A 33 A 12 A 33 A 21 A 33 A 22 A 33 

-, U21^- --, U22 


A13A31 


A13A32 


A 23 A 31 


A 23 A 32 


Let f{xi) - f{xi ± 1) for i = 1,2 be the discrete shift operators in the variables xi 


and X 2 . Introduce the two intertwining operators defined as 


Af = ^ xiT- + ^ X 2 T;^ - ^ (xi + X 2 + p) D, 
Ai 3 A23 ^ A33 




(2.5) 


where D stands for the identity operator. On the bivariate polynomials 
these operators have the actions 


Af Mi^^;,\\xi,X 2 ) = V?L^^M;f^^,„,+i(xi,X 2 ). 


( 2 . 6 ) 


Introduce also the two intertwining operators A(!^ defined in the following way: 

A'J) = Ai^A i 3 + A2iA23 - (Ai^Ai3 + A2^A 23 )n, i = 1 , 2 . (2.7) 

These operators act as follows on the bivariate Meixner pol 5 momials 


AL«Mf„Axi,X 2 ) = 




— 7 (xi X 9 ) 


( 2 . 8 ) 


AL2>Mf,„,(xi,X2) = JjM[^^f^_,(xi,X2). 


The intertwining operators (2.5) and (2.7) can be combined to produce the two 
commuting difference operators Yi and F 2 that are diagonalized by the bivariate 
Meixner polynomials. These operators are defined as 


Yi^AfA^i\ i = l,2. 
Explicitly, one finds 


(2.9) 




AiiA2iA23 


A 


13 


•^ 1 ^X 1 ^*2 + 


AijA2iAi3 


A2iA3fA33^ 


X 2 TX2 


V A23 

''AijA3jAi3) 




AliA3iA33 

Ai3 


xiT, 


XI 


A23 / 

'' A2iA3iA23 
^ A33 


(xi + X 2 + i6) T 


i A33 j 

(xi + X2 + P) + [ A^j xi + A^j X2 + Ag^ (xi + X2 + P)] D. (2.10) 


The eigenvalue equations read 

YiMfln2(^l,X2) = niMl^ln2ixi,X2), T2Mf „2(-^i>-^2) = n2Ml^ln2(xi,X2), (2.11) 

where ni, n 2 are non-negative integers. Let us note that the operators Yi and Y 2 fully 
characterize the polynomials M^^^„^J,x\,X 2 )- 
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3. A discrete and superintegrable Hamiltonian 

We now consider the Hamiltonian obtained by taking the sum of the operators Y\ and 
Y 2 involved in the eigenvalue equations (2.11). We hence define 


Ye^Y^ + Y2. 


(3.1) 


In principle the Hamiltonian (3.1) involves four independent parameters including /3, 
as any matrix A G 0(2,1) depends on three independent parameters. However, it turns 
out that Y6’ essentially depends only on three parameters, including fi. This can be seen 
explicitly as follows. Consider the following parametrization of A in terms of the “Euler 
angles” iy, ^ and ((>: 


Aiy/,^,(p) 


cosh^ 

0 

^sinh^ 


0 

1 

0 


sinh^^ 

0 

cosh^^ 


0 

0 


0 

coshi// 

sinhiy 


0 ' 
sinhiy 
cosiy ^ 


COS0 

-sincp 

I 0 


sincj) 0^ 
cos(f) 0 

0 ij 


(3.2) 


Upon using the expressions (2.10), it is seen that in the parametrization (3.2) the 
Hamiltonian (3.1) has the expression 


JY(y/, (p) - sinh^ if/xiT + cosm if/ sinh^ ^ X 2 T^, 


+ ji- 

X2 ' r —— > Xi-^ X2 

- cosh^ y/ sinh^ ^ (xi + X 2 + j6) - sinh^ yr (xi + X 2 +15) 


n 9 r _ _ 1 9 

cosh ^cosh [xiT^^+X 2 T'^ 2 J + cosh2^cosh y/xi 


+ (cosh^ y/ + sinh^ ^ + cosh^ £, sinh^ y/)x 2 + f5 (sinh^ £, + cosh^ £, sinh^ yr) 


hence the parameter (p does not explicitly appear in The fact that the Hamiltonian 
(3.1) can be presented in terms three independent parameters is a manifestation of its 
superintegrability. Let Jx, Jy and Jz be the operators defined as follows. 

* = iU™A® + AfA<«), 

Jy = i (A"A<_a - Af A'_«), (3.3) 

Jy = i(A™AO>-AfA®). 

The operators Jx, Jy and Jz are constants of motion. Indeed, one can verify by a direct 
calculation that these operators commute with the Hamiltonian (3.1) 


\_Y6', Jx) - 0, [Y6’,Jy) - 0, [^, Jz) - 0. 


The symmetry operators Jx, Jy and Jz satisfy the defining relations of the 5u(2) 
algebra. One has 


\.Jx,Jy) - iJz 


\_Jy,Jz) - iJx 


\-Jz,Jx\ - iJy- 


(3.4) 
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In the realization (3.3), the 5u(2) Casimir operator is related to the Hamiltonian (3.1) 
through 

Jx + Jy + 2 "^ 

The realization (3.3) of the 5u(2) algebra (3.4) and the expression (3.5) of the Casimir 
operator in terms of the Hamiltonian is very close to the Schwinger realization of 
su(2) that one finds when considering the standard two-dimensional quantum harmonic 
oscillator. The SU{2) symmetry (3.4) of the Hamiltonian (3.1) explains why J6’ depends 
on three parameters instead of four: the (p parameter in (3.1) has heen “rotated away” 
from the Hamiltonian hy the choice of parametrization (3.2). 

By construction, the eigenfunctions of the Hamiltonian (3.1) are expressed in 
terms of the two-variahle Meixner polynomials These eigenfunctions 

'F^^^(xi,X 2 ) are labeled by the two non-negative integers N and n and read 

¥^lixi,X2)^Ml^^X-Jxi,X2), 
where n-0,l,...,N and where N - 0,1,2,.... One has 

.^¥^]^(xi,X 2 )^n¥^]^(xi,X 2 ), Jz^fjxyX 2 ) = (n -N/2)wf jxi,X 2 ). (3.6) 

Hence the eigenvalues of are the non-negative integers N -0,1,2,... and are (N+l)- 
times degenerate. The states ^N,n(xi,X 2 ) support (N + l)-dimensional irreducible 
representations of 5u(2). Upon introducing the generators 

J± - Jx ± iJy, 



it is seen from (2.6) and (2.8) that these operators have the actions 


J+'T^;^(xi,X 2 ) 

J-^f,n(xi,X2) 


^n+l)iN-n)¥^]^^^ixi,X2), 

UiN-n + l)wf^^_^ixi,X2). 


(3.7) 

(3.8) 


It thus follows that the two-variable Meixner polynomials Mll^^,.i^ixi,X 2 ) support (K + 1)- 
dimensional unitary representations of 5u(2) where K - ni + n 2 - 

In view of (2.2) wavefunctions 'I'Ar,ra(xi,X 2 ) are not normalized on the infinite grid 
(xi,X 2 ) G N X N with respect to the standard uniform measure of quantum mechanics. 
Properly normalized wavefunctions Y^^^(xi,X 2 ) are obtained by taking 


yfyxi,x2) = 




(xi,X 2 ), 


where (o(xi,X 2 ) is given by (2.3). One then has the orthogonality and completeness 
relations [9] 


^ ^ Y^^^(xi,X2)Yjy, ,j,(xi,X2) = 

X1=0X2=0 


oo N 


L L rZ (Xi,X2)Y^^^(Xi,X 2) = d^^ 

N=0n=0 
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The actions (3.6) and (3.7) on the non-normalized wavefunctions 'I'^^^(xi,X 2 ) can 
he recovered on the normalized wavefunctions Y^^^(xi,X 2 ) hy applying the gauge 
transformation © ej^^^(xi,X 2 ) ^a»“^'^^(xi,X 2 ), where © is either ^ or any one of the 

symmetries Jx, Jy, Jz- 

Below are illustrated some of the wavefunctions amplitude I Y^^^(xi,X2)| for various 
values of the parameters i/a, 6 and The model is defined on the N x N grid hut only 
the grid 75 x 75 is shown. It is seen that the particle is localized close to the origin. The 
energy level N prescribes the number of “nodes” in the wavefunction amplitudes and 
the parameter (5 is related to the spatial spreading of the amplitude. 





Figure 1: Wavefunction amplitudes for ^ = 0.8, y/ - 0.8, (p '■ 


|Y 


(/ 5 )| 

0 , 0 ' 


middle line: 


top line: IY^^q 


n/4 and p - 15. Bottom line: 



Figure 2: Wavefunction amplitudes for ^ = 0.8, y/ - 0.8, (p - nl4 and /3 = 28. Bottom line: 
lYjfJl, middle line: lY^^Jl, lY^^Jl, top line: |Y^^J|, |Yg|, |Y^^^|. 
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Figure 3: Wavefunction amplitudes for ^ = 0.5, f - 0.8, (p - n/A and /3 = 15. Bottom line: 


Imiddle line: top line: ||Y 


2 , 2 ' 



Figure 4: Wavefunction amplitudes for ^ = 0.8, p/ = 0.8, (p ■ 


|Y 


0 , 0 ' 


middle line: IY^^^qI, 


top line: IY^^q 


0 and 13-15. Bottom line: 


4. Continuum limit to the standard oscillator 

It will now be shown that in the continuum limit, the model governed by the 
Hamiltonian (3.1) tends to the standard isotropic quantum oscillator in two dimensions. 
The continuum limit from the two-variable Meixner polynomials to a product of two 
Hermite polynomials will be considered first. The explicit limit of the Hamiltonian 
(3.1) and of the constants of motion (3.3) will then be investigated. 
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4.1. Continuum limit of the two-variable Meixner polynomials 


Consider the generating function (2.1) of the two-variahle Meixner polynomials (2.4) in 
the parametrization (3.2). Upon writing 




a 




21 


2 l 


22 


22 


V?’ ^ V?’ 

in left-hand side of ( 2 . 1 ) and taking the limit as f ^oo using the standard result 


(4.1) 


lim 

k^oo 




one finds that 


lim 

)3^oo 


^ ^ A 31 Zl ^ A 32 22 


-X 1 -X 2-/3 


A33 A33 ^ J 

( All 21 Ai2 22 \ 

^ Ai3^ 

- exp [- 2 i(acos(^- 6sin0)] 

X exp [ - 22(0 sin (p + b cos <p)] 


xi / . . \X2 

, L A 21 21 A 22 22 

1 + ^-+ ^-;= 


A23^ A23^) 


1 + 


cos<p sin 0 


a 


-21 + 


a 


-22 


XI 


sin (p cos (p 


b 


-21 + 


b 


-22 


X2 


Upon defining 


C'„i,„2(-^i,-^ 2)= lim (xi,X2), 

p—*oo 

and taking the limit as 00 with (4.1) in the right-hand side of (2.1), one finds that 


exp [- 2 i(acos(^- fesin^))] exp [- 22(0 sin (^-l- &cosd>)] 


'' cos(p sin(p 


1 + 


a 


-21 + 


a 


-22 


'' sin 0 cos<p 00 00 


b 


21+ —^22 = L L Cm^^n2(xi,X2) 


?ll= 0«2 = 0 


~ni n2 

\/ni\n2\ 


. (4.2) 


It is seen that the pol 5 momials correspond to the two-variahle Charlier 

polynomials [ 8 ]. If one uses the parametrization (4.1) and takes the limit /3 ^ 00 in the 
weight function (2.3), one finds the two-variable Poisson distribution 


hm w(Xl,X 2 ) = e ^ ^-;-;—, 

p^oo Xi!X2! 

and the orthogonality relation ( 2 . 2 ) becomes 


(4.3) 


E I 

X1=0X2=0 


xi!x 2 ! 


^n\,n2 (xi,X2)C,j'^ ,j^(xi,X2) ^ni,n'^^n2,n.2' 
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It is directly seen from (4.2) and the standard generating function for the one-variable 
Charlier polynomials [18] that when (p-0, one has 




(f>=0 


(_ 

anibn2 


Cniixr,a^)Cn^{X2-,b‘^), 


(4.4) 


where C„(x;a) is the one-variable Charlier polynomials. Thus, using the standard limit 
from the one-variable Charlier polynomials to the one-variable Hermite polynomials, 
one can set 


xi-V2axi + a‘^, X2 = n/26x2 + 6 ^, <p-0, 


(4.5) 


in (4.2) and take the limit as a ^ oo and 6 ^ oo to find 


lim lim 

a^oob^oo 


+ V2xizi^--^ + V2 


X2 Z2 


Upon comparing with the well-known generating function for the Hermite polynomials 
[18], one finds that 


lim lim Cni,n 2 ^x\,X 2 ) 

a^oob^oo 


0 = 0 


\j2^^+^‘^ni\n2\HnSxi)Hn.^{X2), 


where H„(x) are the standard Hermite polynomials. With the parametrization (4.5), it 
is easily shown using Stirling’s approximation that the bivariate Poisson distribution 
appearing in (4.4) converges to the normal distribution 


lim 6 "“' 

a^oo 


(o^ 

xi! 


e 


\/n 


In summary, the wavefunctions (3.9) of the discrete two-dimensional system governed 
by the Hamiltonian (3.1) tend to a separated product of two univariate Hermite poly¬ 
nomials in the combined limiting process (4.1) and (4.5). This motivates calling (3.1) a 
discrete oscillator. For other limits of bivariate orthogonal polynomials, see [1]. 


Remark 

It is not needed to take 0 = 0 in the second limiting process involving the two-variable 
Charlier polynomials. If one keeps (p arbitrary and performs the change of variable 
(4.5) and takes the limit as a ^ oo and b oo, one simply finds a product of Hermite 
polynomials in the rotated coordinates xi = cos 0 xi - sin 0 X 2 and X 2 = sin 0 xi -l- cos 0 X 2 . 


4.2. Continuum limit of the raising j lowering operators 

Let us now examine the combined limiting procedures of the preceding Subsection and 
its effect on the defining operators of the discrete oscillator model defined by (3.1). We 
first consider the raising and lowering operators (2.5) and (2.7) of the bivariate Meixner 
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polynomials. Under the gauge transformation = u»^^^(xi,X 2 ) A^^ o ^^^(xi,X 2 ), these 
operators have the expressions 

A + ^ = Aii + A2i \/xi +X2 + /3 D, 

A^i^ = All \/^i + 1 + A2i \/x2 + 1 — A^i \/xi + X2 + P D, 

for i = 1,2. On the wavefunctions (3.9), one has 

A<« Y^;'\xi,X 2 ) = 

Af Y^;'\xi,X 2 ) = ^N-n + 1 A%\^^^{xi,X 2 ), 

and 


A<_« Y^J„(xi,X 2) = ^ Y;f_\^;^_^(xi,X2), 
A® Y^J„(xi,X 2) = Y^j|^^_\(xi,X2). 


Upon taking the parametrization (4.1) and taking the limit as ^ ^ oo, the raising 
operators become 


ar^’ - lim A'^^^ = cos (p T - sin (p sjx 2 T - (a cos (p-b sin0) [ 

P^CX) ^ ^ 


a^^’ - lim A''^ - sincp T + cos cp T - (a sin0 + b cos cp) I 

P^OO ^ ^ 


and the lowering operators become 


a 


_ = lim A_^ = cos(p y/xi + 1 - sin0 \/x 2 + l - (acoscp- bsin(p) 


X2 


a 


P^OO 

- lim A^^^ = sin0 y/xi + 1 + cos0 \/x 2 + l - (asin0 + bcoscp) I 


P^OO 


(4.6) 


(4.7) 


A direct calculation shows that these operators satisfy the commutation relations 

[a^l’,af] = 6ij, = 0 , , a^'^] = 0 . 

Upon setting xi = \/2axi + a^ and X 2 = V2bx2 + b^ as in (4.5) and taking the limit as 
a ^ oo and 6 ^ oo, it is easily seen that the operators (4.6) and (4.7) become the rotated 
creation/annihilation operators 

^ cos 0 a^ -sin^Ug, ^ sin 0 a^ + cos 0 a 2 , 
a^}^ cos0ai -sin0a2, a^3^ sin0ai + cos0a2. 


where 


ai = 


Xi + dxi 

~vr 



Xi ^Xi 


are the usual creation/annihilation operators. It immediately follows that 
in the continuum limit described above in a two-step process, the gauge- 
transformed Hamiltonian (3.1) and constants of motion (3.3) obtained through 
a»(xi,X 2 )^^^^e>“^^^(xi,X 2 ) tend to the standard two-dimensional oscillator Hamiltonian 
and the 5u(2) generators in the Schwinger realization. 
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5. Conclusion 

In this paper, we have introduced and described a discrete model of the oscillator in 
two-dimensions based on the bivariate Meixner pol 5 momials. We have shown that this 
system is superintegrable and that it has the same symmetry algebra as its continuum 
limit, the standard isotropic oscillator in two dimensions. We have established that the 
two-variable Meixner polynomials form bases for irreducible representations of 5u(2). 
We have detailed the limiting processes by which the two-variable Meixner pol 5 momials 
tend to the bivariate Charlier polynomials and by which the latter tend to a product of 
standard Hermite polynomials. 

In the present paper we have considered for simplicity the two-dimensional case. 
However, since the theory of multi-variable Meixner is now well established, it is 
clear that the model can be generalized to any dimensions to give a d-dimensional 
model of the harmonic oscillator with the same 5u((i) symmetry as the standard 
quantum oscillator in d-dimensions. Another possible generalization would be to 
consider, instead of (3.1), a discrete anisotropic oscillator with a Hamiltonian of the 
form -Yx + aY 2 - It is clear that for rational values of a this model would still be 
superintegrable, but would exhibit a higher order symmetry algebra. The case a = 2, 
which would correspond to a discrete 2 :1 anisotropic oscillator, is of particular interest 
in that regard [7]. 
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